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The stability of solutions to evolution equations with respect to small stochastic perturbations 
is considered. The stability of a stochastic dynamical system is characterized by the local stability 
index. The limit of this index with respect to infinite time describes the asymptotic stability of a 
stochastic dynamical system. Another limit of the stability index is given by the vanishing intensity 
of stochastic perturbations. A dynamical system is stochastically unstable when these two limits 
do not commute with each other. Several examples illustrate the thesis that there always exist 
such stochastic perturbations which render a given dynamical system stochastically unstable. The 
stochastic instability of quasi-isolated systems is responsible for the irreversibility of time arrow. 
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I. INTRODUCTION 

o 

' Evolutional processes of nature are described by differential equations that, in general, are equations in partial 
derivatives. A set of such partial differential equations constitutes an infinite-dimensional dynamical system. Under 
a physical system one implies en ensemble of objects whose behaviour is of interest. The evolution of a given physical 
system is characterized by the related dynamical system. Among physical systems, one distinguishes isolated systems 
as opposed to open systems. The evolution of the isolated physical systems is governed by deterministic laws, that 
is, by deterministic equations, not containing random variables. While open physical systems, generally, deal with 
stochastic equations, where random terms represent the interaction with surrounding. 

Solutions to differential equations can be either stable or unstable. There are methods for analyzing the stability of 
, solutions for a given dynamical system, either deterministic [1-3] or stochastic [4]. Here we address another problem, 
that of stability of a deterministic dynamical system with respect to small stochastic perturbations. This problem is 
not only interesting by itself but it is of fundamental importance with regard to the question: How adequately the 
■ notion of isolated systems represents the physical reality? 

As is evident, the notion of an isolated system is an abstraction. In fact, no real system can be completely 
isolated from its surrounding. This point has been repeatedly emphasized in literature [5-9]. And the impossibility 
of ideally isolating macroscopic systems from their environment is considered as being intimately related with the 
irreversibility of time [10,11]. Moreover, it has been stressed [12,13] that the concept of an isolated system is logically 
self-contradictory by its own. This is because to realize the isolation, one has to employ technical devices acting 
on the system; and to ensure that the latter is kept isolated, one must apply measuring instruments perturbing the 
system. The preparation and registration processes disturb the system dynamics [14]. In this way, there exists an 
accepted understanding that any considered physical system is never absolutely isolated but is subject to, probably, 
weak but, generally, uncontrollable random influence from the environment. Even if this influence is quite weak, its 
very existence is of principal importance, for explaining the irreversibility of time arrow. 

It is worth noting that the irreversible behaviour of macroscopic systems is often attributed to internal chaotic nature 
of microscopic dynamics (see discussion in [15]). However, not all physical systems display chaotic behaviour. Many of 
them are perfectly governed by rather simple deterministic laws, with no signs of chaos. Nevertheless, the time arrow 
is well defined for any system, including very simple and not chaotic ones. What is more, the recent developments in 
dynamical theory, as reviewed by Zaslavsky [16], show that chaotic dynamics in real systems does not provide finite 
relaxation time to equilibrium or fast decay of fluctuations, and that chaotic systems are not completely random in 
the sense originally postulated for statistical systems. Therefore the presence of random environment, though very 
weak, seems to be crucially important for interpreting fundamental notions in the behaviour of real physical systems. 
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From another side, there is a common belief, based on practical experience, that physical systems can, with a very 
good accuracy, be isolated and can be described by deterministic equations, while the random influence of surrounding 
may be neglected. Thus, there exists an apparent contradiction between the principal necessity of allowing for random 
perturbations influencing any real system and the practical possibility of neglecting such perturbations, treating a 
system as isolated. 

This contradiction is resolved in the present paper by putting the problem on a firm mathematical footing. The 
concept of quasi-isolated systems is defined. It is shown that such systems, generally, are unstable with respect to 
infinitesimally small stochastic perturbations. At the same time, for a finite temporal period, these systems can be 
treated as approximately isolated. 



II. STABILITY OF STOCHASTIC SYSTEMS 

Let a continuous variable x € D denote a set of spatial coordinates pertaining to a domain D and let t € R+ denote 
time. Suppose a stochastic field £(t) is defined. In general, the latter is a set of stochastic functions £i(x,t), with 
z = l,2,.... Throughout the paper, we shall use the matrix notation [17] making it possible to express the following 
equations in a compact form. Thus, the stochastic field £(t) = [£»(ar, t)] is considered as a column with respect to both 
i = 1, 2, ... as well as x € D. The dynamical state ?/(£, t) = [yi(x, £, t)j is also a column with respect to i and x, as is 
the velocity field v(y, £, t) = [vi(x, y, £, t)]. The set of evolution equations, defining a dynamical system, in the matrix 
notation reads 

j t y(ti,t) = v(y,ti,t). (1) 



This is complimented by an initial condition 



y(£,0) = y(0), (2) 



implying the set 

yi(x,Z,0) = yi (x,0) (i = l,2,...) 

of the related initial conditions. The averaging over the stochastic field £(f) is denoted by the double angle brackets 
as 

y(t) =< y(U) » , (3) 

which assumes the family of the functions 

yi(x,t) = 4Zyi(x,£,t)», (4) 

with i = 1,2,.... 

In the stochastic equation (1), the velocity field v(y,£,t) may, in general, contain differential as well as integral 
operations. To solve Eq. (1) means to find the averaged solution (3). Stochastic differential equations, as is known 
[18], can be defined either in the sense of Ito or in the sense of Stratonovich. In what follows, the latter definition 
will be employed, which permits simpler calculations and is better motivated physically [19]. It is also possible to use 
the stochastic expansion technique [20,21], presenting the stochastic field as an expansion over smooth functions of 
spatial and temporal variables with random coefficients. This method enables the usage of the standard differential 
and integration analysis. The final results of the expansion technique coincide with the corresponding expressions 
obtained by means of the Stratonovich method. 

The local stability of a dynamical system can be characterized by the local stability index 

^InsupM (5) 

6y(0) \m0)\ 

which describes the maximal deviation of the averaged trajectory at time t after an infinitesimal variation of the 
initial conditions. Such a deviation, according to Eq. (5), corresponds to the law 

\Sy(t)\ ~ \Sy(0)\e°W , ( 6 ) 
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from where it is evident why cr(t) is called the stability index, or stability exponent. From this definition, one can 
immediately conclude that the admissible local properties of motion are classified as: 



a(t) < (locally stable) , 

a(t) = (locally neutral) , 
a(t) > (locally unstable) . (7) 
The asymptotic Lyapunov stability corresponds to the terminology: 

lim a(t) = — oo (Lyapunov stable) , 

t — >oo 

lim a(t) > — oo (Lyapunov unstable) . (8) 

t — >oo 

And in the language of the Lagrange stability of motion, one has: 

sup a(t) < oo (Lagrange stable) , 
t 

supcr(t) = oo (Lagrange unstable) . (9) 
t 

The limit 

A = lim - a(t) (10) 

t — >oo t 

corresponds to the largest Lyapunov exponent. One tells that the motion is asymptotically stable if A < 0, neutral 
when A = 0, and unstable if A > 0. 

The usage of a local characteristic of motion, such as the local stability index (5), provides us an essentially richer 
information on temporal dynamics than the largest Lyapunov exponent (10) defined for the limit t — > oo. First of 
all, this is because many dynamical systems possess a rather complicated structure of their phase space resembling a 
topological zoo, consisting of domains of chaotic dynamics as well as of regions of regular motion, containing manifolds 
of wandering trajectories as well as trapping islands. As a result of this, the fine local properties of orbits play a 
leading role, while such a fairly rough characteristic as the limiting Lyapunov exponent is less important [16,22]. 

Moreover, the asymptotic divergence of trajectories of stochastic dynamical systems is not compulsory exponential 
[4], because of which making use of only the limiting Lyapunov exponent (10) may result in the loss of information. 
For example, the divergence of trajectories can be of power law 

\Sy(t)\ ~ |5j,(0)|t" . 

Such power laws are typical for weakly disordered systems [23] exhibiting mid-range order [24]. In that case, the 
local stability index (5) behaves as a(t) ~ (3 hit, which can be either positive or negative depending on the sign of [3. 
Respectively, the motion is either stable or unstable. While, according to the Lyapunov exponent (10), which is A = 0, 
the motion is neutral. Another example has to do with the divergence of trajectories by the stretched exponential law 

\6y(t)\ ~ \5y{0)\exp(KtP) , 

with < j3 < 1, which is also quite ubiquitous in disordered systems. Then the local stability index (5) is a(t) ~ nt^ , 
which again can be either positive or negative depending on the sign of k, hence, the motion is either stable or 
unstable. And the limit (10) is again zero, classifying the motion as neutral. 

Instead of the asymptotic Lyapunov exponent (10), one could define the local Lyapunov exponent [25,26] as 

A(t) = \ a(t) . 

However, for what follows, the usage of the local stability index (5) is more convenient. 

One more advantage of employing a local characteristic of stability is that the limit (10) for many complex systems 
is technically unachievable. Then the local index (5) is the sole available quantity that can be actually calculated. 
Such a situation is typical for complicated nonlinear equations that can be treated only numerically [27], for the 
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analysis of various time series that are always finite [28] , and for the dynamical representation of perturbation theory 
where it is practically feasible to calculate only a finite number of terms [29-32]. 

The local stability exponent (5) can be expressed through the multiplier matrix M(t) = [M i} -,{x, x' , t)\ with the 
elements 



Mij(x,x',t)= ^ Y„/ rVi • (11) 



From this definition, it follows that 

M ij (x,x',0)=6 ij 6(x-x') , (12) 

where Sij is the Kroneker delta and S(x) is the Dirac delta-function. Writing the variation of the averaged dynamic 
state as 

5y(t) = M(t)5y(0) , (13) 

we see that 

sy(o) \m0)\ 

with the spectral norm of M{t) being assumed. Therefore the local stability exponent (5) is 

a{t)=\n\\M{t)\\. (15) 
Thus, to analyse the stability of motion, we need to know the multiplier matrix (11). 

III. STOCHASTIC MULTIPLIER MATRIX 

What we are actually given is the stochastic equation (1) defining the stochastic dynamic state y(£,t), whose 
variation 

5y(£,t) = M(£,t)5y(0) (16) 
over the initial conditions involves the stochastic multiplier matrix M(£,t) — [Mij (x, x' , £, t)] with the elements 

^•'•m-Wm- (17) 

For the latter, one has the initial condition 

M ij (x,x',£,0) = 5 ij 5{x-x') . (18) 

The multiplier matrix (17) is connected with the stochastic Jacobian matrix J(£,t) = [Jij(x,x',£,t)] with the 
elements 

Mx^i) 

j lA x,x,^,t)- Syj{x ,^ jt) ■ U»J 

The variational differentiation of Eq. (1) gives the equation 

j t M(U) = J(U)M(Z,t) (20) 

for the multiplier matrix (17). The initial condition for this equation is Eq. (18). 

Since the evolution equation (1) represents a set of partial differential equations, one has to define as well boundary 
conditions. The latter can be written in the general form 

b(y,t,t) = (xedB), (21) 
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where SO is the boundary manifold of the domain © and b(y, £, t) — [bi(x, y, £, t)] is a boundary vector. Defining the 
boundary matrix B(£,t) — [Bij(x, x' , £, t)] with the elements 



B-ij (x, x , i) 



<>k(x,y,£,t) 



(22) 



5yj(x',£,t) 

and accomplishing the variation of Eq. (21), we get the boundary condition 

B(£, t)M(£, i) = (x e dB) (23) 

for the multiplier matrix. 

As an illustration, we may offer the often met form of the boundary conditions 

(i + C^)yi(*,£,*) = /i(t) (xedB), 

where ( is a parameter and fi(t) is a given function. The variation of this condition results in the equation 

(l + {-j£)M ij (x,x',Z,t) = (xedB), 

demonstrating a particular case of the boundary condition (23). 

For the multiplier and Jacobian matrices, one may employ different representations. To this end, let a set {(fi n (t)} 
of the columns ip n (t) — [<p n i(x,t)] be given, forming an orthonormalized complete basis, 

<Pm(t) <Pn(t) = 5 mn , ^ ^(t) = * ; 

n 

where 1 = [SijS(x — x')] is the unity matrix and n is a labelling multi-index. To pass from the x-representation to 
n-representation, we define 



M mn (£,t) = <P+(t) M &t)<Pn(t) , Jmn&t) = >fi+(t) J t)<p n {t) . 



(24) 



Recall that the matrix notation [17] is used here, according to which, for instance, the action of the multiplier matrix 
on <p n (t) is the column 



M(£,t)<p n {t) 



f Mij(x,x',€,t)<p n j(x',t) dx' 



Equation (20) for the multiplier matrix in the new representation reads 



where the relation 



Jmk{^t)M kn {^t) + M mk {^t) V +{t) - <p+(t) M kn (Z,t) 



d<p n (t) +M dif k (t) 



(25) 



-^— iPn{t ) + iPm (t)^ r =0, 
following from the normalization condition, is used. And from Eq. (18), we have the initial condition 

for Eq. (25). The multiplier matrix enjoys several useful properties. 

Proposition 1. If the dynamical state y(£,t) can be presented as an expansion 

2/(£,i) = ^>„(£,i) ip n (t) + f(t) , 



(26) 



(27) 
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over a basis {f n (t)} and f(t) — [fi(x,t)] is a column of functions not depending on the initial state y(0), then the 
multiplier matrix has the form 

M(£, t) = J2 M£> *) fn(t) <p+(0) , (28) 

n 

in which 

Proof. The variation of the expansion (27) gives 

$yt(x,Lt) = \p 5c n (£,t) 6c n (£,o) 

5 yj (x',0) ^6c n (£,0) 6 yj (x',0) ^ MlX ' >■ 
At the same time, from Eq. (27) we have 

c n (Z,t)=<pt(t)y(Z,t)-<p+(t)f(t) . 

From the latter equation, we get 

5c n (lj,0) , 
%(a/,0) ^ (X ' Uj - 

Using this and invoking the definition (17), we obtain the form (28) with notation (29). 

Remarks. Although the basis {<fn(i)} is assumed to be orthonormalized, but the vectors (p m (ti) and (p n (t2) at 
different times t\ ^ t 2 arc not necessarily orthogonal, so that, in general, 

Neither ip n (t) nor ip n (0) are necessarily the eigenvectors of the multiplier matrix, for which we have 

M(£,t) (f n (0) = H n {S,t)<p n (t) . 

Only when (p n (t) = ip n does not depend on time, then tp n is an eigenvector of M(£,i) and /z„(£,i) is its eigenvalue. 

Proposition 2. Suppose the multiplier matrix M(£,t) possesses eigenvectors (p n (t) forming a complete orthonor- 
malized basis. Then the related eigenvalues, given by the eigenproblem 

M(£,t)(p n (t) = M& *)¥>»(*) , (30) 

can be presented as 

/x„(£,i) =exp( / Jnnitt') dt'X . (31) 



Proof. With (p n (t) being the eigenvectors of the multiplier matrix, the elements of the latter, defined in Eq. (24), 
are 

(£,t). (32) 

Substituting this into Eq. (25) yields 

Smn J t Hn(Lt) = Jmn(U)»n(U) + [M£> *) ~ M£> . (33) 



When to = n, the latter equation gives 



^M£,t) = -M£,*K(£,t), ( 34 ) 
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while for m ^ n, it results in 



Solving Eq. (34), with the initial condition 



, dip n (t) 

v m («) ~ ir 



Mf.0) = l. (35) 



we come to the eigenvalue (31). 

Remarks. From the eigcnproblem (30), one gets the representation 

M(£,t) = J2fi n (Z,t)<p n (t)<p+(t) (36) 

n 

for the multiplier matrix. The eigenvectors of the latter are not necessarily the eigenvectors of the Jacobian matrix 
(19). Hence the form J mn (£,,t), defined in Eq. (24), is, in general, nondiagonal. 

Proposition 3. Assume that a complete orthonormalized basis {<£«(£)} is such that 

^ - (m + n) . (37) 

Then ip n (t) are the eigenvectors of the multiplier matrix M(£,t) if and only if they are also the eigenvectors of the 
Jacobian matrix J(£,i). 

Proof. Let condition (37) hold. Then Eq. (25) becomes 



^ M rnn (Z,t)=J2 J rnk(Z,t)M kn (Z,t)+M mn (Z,t) 



dt 

k 



■ . d<p n (t) . dy m (t) 



(38) 



If (/5„(i) are the eigenvectors of M(£,t), that is, the form (32) takes place, then Eq. (38) reduces to 

3mn l^ni^jt) Jmn (£; fyftn {£,■> t) , 

from where it is clear that 

Jmni^it) — fimnJnni^it) ■ (39) 

Hence, <p n (t) are the eigenvectors of J(£, t). 

Conversely, if (p n (t) are the eigenvectors of J(£,t), so that Eq. (39) holds true, then solving Eq. (38) yields 



dt'\ . 



M mn %t) = M„ m (e,o)exp|^ \j mm {a,t')+vt{t') ^(f) d ^f 1 

In view of the initial condition (26), this results in 

M mn {tt) = <5 mn exp J„ n {Z,t') di'J , (40) 

which tells us that (p n {t) are the eigenvectors of M(£,t). 

Remarks. As follows from Eq. (40), the eigenvalues of the multiplier matrix are given by expression (31). A 
simple example, when condition (37) is valid, is the case of a stationary basis {</?«}, with ip n (t) — ip n not depending 
on time. 

Comparing Eqs. (3), (13), and (16), we see that 

M{t) = «M(£,i) > . (41) 
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Therefore, if M(£,t) possesses eigenvectors (p n (t), then the matrix (41) satisfies the eigenproblem 

M(t)<p n (t) = Hn(t)<p n (t) (42) 

with the same eigenvectors and the eigenvalues 

M*) = «M&*)»> (43) 

which have the property 

MO) = 1 • (44) 

With the spectral norm 

\\M(t)\\ =sup|Mt)| . 

n 

the local stability exponent (15) becomes 

a(t) = In sup | < Hn(£,t) > | • (45) 

n 

In this way, the problem of analyzing the stability of a stochastic dynamical system is connected with finding the 
eigenvalues of the stochastic multiplier matrix. 

IV. CONCEPT OF QUASI-ISOLATED SYSTEMS 

As is discussed in the Introduction, no real physical system can be completely isolated from its surrounding. The 
latter can be modelled by stochastic perturbations of the system dynamics. To stress that the amplitude of the 
stochastic perturbation is small, it is convenient to include explicitly a small factor a in front of the stochastic field 
£(i). So, instead of Eq. (1), we shall write 

(46) 

The factor a = oi\ + ia 2 is assumed to be complex, with its real part a\ = Re a and imaginary part a 2 = Im a. The 
complex value of the factor a makes it possible to simulate random fluctuations of different physical quantities, such 
as energy and attenuation or density and phase. If a = 0, there are no stochastic perturbations, and one returns to 
a deterministic dynamical system. When stochastic fields are switched on by means of a ^ 0, we have a stochastic 
dynamical system, whose local stability is characterized by the stability exponent (45) that takes the form 

<r(a,t) = In sup | < M^f.*) » I > ( 47 ) 

n 

where the dependence on the switching factor a is explicitly shown. 

The stability exponent (47) describes the stability of a stochastic dynamical system with respect to the infinitesimal 
variation of initial conditions. For correctly defining the notion of a quasi-isolated system, it is also necessary to 
consider the stability with respect to infinitesimal stochastic perturbations. This implies that, after analyzing the 
stability of the stochastic system by means of the stability exponent (47), we should set a — > 0. Since a is complex- 
valued, the limit a — > means that both its real and imaginary parts tend to zero: a.\ — * and a 2 — > 0. Among 
all admissible ways of tending to zero for a — > 0, it is necessary to chose that one providing the maximal value for 
the exponent (47), in agreement with its definition as characterizing the largest deviation of the trajectory. The so 
defined limit a — > will be denoted as 

lim a(a, t) = sup lim <r(a,t) . (48) 

a^O a \ a \— >0 

In the stability analysis with the help of the stability exponent (47), an important part is the consideration of the 
asymptotic stability, when t — > oo. This limit may, in general, not commute with the limit a — > 0. Therefore, an 
important step is to study the commutativity of these limits, characterized by the commutator 

[ lim , lim ] = lim lim — lim lim . 

a— >0 t-^oo a^0 t^oo t—>oo a^Q 
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The content of this section can be summarized by formulating the following definitions. 

Definition 1. A physical system is called quasi-isolated if its evolution is described by the stochastic dynamical 
system (46) with infmitcsimally small stochastic perturbations. 

Definition 2. A quasi-isolated system is stochastically stable when 

[lim, lim ]a(a,t) = . (49) 

a^O t— >oo 

Definition 3. A quasi-isolated system is stochastically unstable if 

[lim, lim]a(a,t) ^ . (50) 

Note that the inclusion of stochastic fields in the evolution equations can be realized in different ways. Hence, 
in principle, one could consider the stochastic stability with respect to each of particular ways. A quasi-isolated 
system may turn to be stochastically stable with respect to some of perturbations but unstable with respect to others. 
However, the kind of action of random environment on a quasi-isolated system is, by assumption, unpredictable. 
Therefore, it is not sufficient to limit ourselves by only some ways of including stochastic perturbations, which would 
result in the analysis of partial stochastic stability. But, in order to make conclusion on the general stochastic 
stability of a quasi-isolated system, one must analyze all qualitatively different admissible ways of including stochastic 
terms in the evolution equations. Fortunately, there are just two main qualitatively different types of random noise, 
multiplicative and additive. 

In the following sections, the preceding ideas will be illustrated by concrete examples. Since the perturbing influence 
of surrounding may be caused by many independent random sources, their action, according to the central limit 
theorem, can be modelled by the Gaussian white noise [18]. For the convenience of the reader, the basic properties of 
this noise, which will be repeatedly used throughout the paper, are listed in short in the Appendix. 



V. IMPORTANCE OF MULTIPLICATIVE NOISE 



One may notice that additive noise cannot lead to stochastic instability. Really, let the velocity field in Eq. (1) be 
a sum v(y, £, t) = v\ (y, t) + V2{£, t) of two terms, the first of which does not depend on the stochastic field £(i), while 
the second does not include the dynamic state y. Then the Jacobian matrix (19) is defined only through the variation 
of vi and does not depend on v 2 - Therefore the solution of Eq. (20) for the multiplier matrix also is independent 
from v 2 , which means that v 2 does not influence the properties of the multiplier matrix, hence, does not change the 
type of stability. 

But the multiplicative noise can strongly influence the stability property To illustrate this, let us consider the 
evolution equation (46) with the velocity field 

v(y,aZ,t) = f(t) + at(t)y(at,t) , 

where f(t) is a given function and £(£) is a Gaussian white-noise variable with the properties described in the Appendix. 
The equation (46), 

§=/(t)+a£(t)y, (51) 

determines the evolution of a one-dimensional dynamical system. In this case, the Jacobian matrix (19) reduces to 
the function 

J(a£,t) = a£(t) . 

According to Eq. (31), this gives the multiplier 

/i(a£,t) = cxp 

The same form (52) could be obtained from the direct variation of the solution 

y«, t) = t/(0) exp ja f £(t') dt'j + f f(t') cxp ja J £(t") dt" J dt' . 
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For the stability index (47), we find 

a(a,t) = (al-aDjt, (53) 

where a\ = Re a and = Im a. Keeping in mind the definition (48), according to which the stability index is to 
be maximized with respect to oti and c*2, under the given modulus \a\ 2 = a\ + we see that sup Q (a 2 — a 2 .) equals 
| a | 2 = a\. Therefore the index (53) can be written as 

a(a,t) = \a\ 2 -yt . (54) 

From here it follows that the limits 

lim lim a(a, t) = , lim lim a(a, t) = oo (55) 

do not commute with each other. This implies that the quasi-isolated system, whose evolution is given by Eq. (51), 
is stochastically unstable. 



(57) 



VI. OSCILLATOR IN STOCHASTIC BACKGROUND 

Many physical processes are presented by oscillatory motion. It is, therefore, illustrative to consider a quasi-isolated 
system described by a harmonic oscillator subject to the action of a weak external noise. Let the evolution equation 
(46) have the form 

^ = iuy + a£(t)y , (56) 

where the oscillator frequency u> is real. Here the real part of a corresponds to the noisy attenuation-generation 
process and the imaginary part of a described the noise of frequency. 
For this one-dimensional case, the Jacobian matrix (19) is the function 

J(a£,t) =itu + a £(t) . 

In view of Eq. (31), the multiplier is 

At(a£,i) = cxpjiwt + a J f(t') di'j . 
The same expression (57) also follows from the variation of the solution 

y(a£,t) =j/(0)cxp jiwt + ajf £(t') ctt'J . 

The stability index (47) is 

a(a,t) = (aj-aD-ft , (58) 

where the properties of the white noise from the Appendix are used. 

If the influence of the random noise is removed before the temporal limit, that is, a — > 0, then, for any choice of a.\ 
and «2) we have 

lim lim a(a, t) = , (59) 

t — >oo a — >0 

which corresponds to the neutral motion. However, the situation is different if the limit t — > oo is taken first. Then, 
maximizing the factor (58), in agreement with definition (48), as is explained in the previous section, we get the form 
(54). As s result, 

lim lim a(a,t) — oo . (60) 

The noncommutativity of the limits (59) and (60) shows that the oscillatory motion is stochastically unstable. 

This means that for a finite time, such that |a| 2 7i <C 1, the system with an oscillatory evolution can approximately 
be treated as isolated. But there always exists such a weak random noise that makes the system unstable for sufficiently 
long times. 
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VII. STOCHASTIC DIFFUSION EQUATION 



Consider the diffusion equation 



| = [* + <*(*)] g, (61) 

in which the diffusion constant D > is subject to weak random fluctuations. For any given finite interval, the spatial 
variable x can always be scaled so that to be defined on the unity interval. Thus, we assume that x € [0, 1]. Equation 
(61) is complimented by the initial condition 

y(x,a£,0)=v(x,0), (62) 
with a given function y(x,0), and by the boundary conditions 

y(OM,t) = bo, y(lM,t) = h, (63) 

where bo and b\ are constant. 

For Eq. (61), the Jacobian matrix (19) is 

d 2 

J(x, x', aO = [D + a£(t)} S(x - x') . (64) 

The boundary conditions (63) lead, according to Eqs. (22) and (23), to the boundary conditions 

M(0, x', a£, t) = M (1, x', a£, t) = (65) 

for the multiplier matrix. 
Solving the eigcnproblcm 

l 

J(x,x',a£)ip n (x') dx' = J n (a(,)ip n (x) (66) 



for the Jacobian matrix (64), with the boundary conditions 

<p n {0) = <p n {l) = , (67) 



we find the eigenvalues 



and the eigenfunctions 



where 



J„K) = -[D + a^t)} k 2 n (68) 
fn{x) = V2 sink n x , (69) 

k n = im (n= 1,2, ...,JV-» oo) . (70) 

The eigenvectors ip n = [ip n (x)], being the columns with the elements (69), are stationary. Hence, they satisfy 
condition (37). Then, by theorem 3, the multiplier matrix possesses the same eigenvectors ip n , with the eigenvalues 
(31), where J nn = J n . Taking account of Eq. (68) yields 

ft 

,2. 2 / t (*)\ ,H' 



M«e,*)=exp|-D^t-a^ J £(t')dt'j . (71) 
Note that the solution to Eq. (61) reads 

oo 

(a£,t)<p n (x) + f(x) , 



n=l 
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where 



c n = / [y(x, 0) - f{x)]ip n (x) dx , f(x) = b Q + (bi - b )a 
Jo 



The form of this solution is that of the expansion (27) in theorem 1, because of which the multiplier matrix could be 
found by means of this theorem. 

Averaging Eq. (71) over the stochastic field (see Appendix), we get 

| < At„(a£, > I = CX P {- Dk lt + oPkfat) , (72) 
where a is real. Hence, the stability index (47) becomes 

a(a,t) = sup (-Dklt + a 2 k 4 n jt) . (73) 

n 

Taking into consideration Eq. (70), this gives 

. .,_( -DttH (a = 0) 
~ \ ct>(*N)*>yt (a^O) 

with N — > oo. 

In this way, we have 

lim lim a(a, t) = — co , (74) 

which means that in the absence of any stochastic perturbations the motion would be stable. However, if infinitcsimally 
small stochastic perturbations are present, then 

lim lim a(a,t) — co , (75) 

and the motion is stochastically unstable. This case serves as a good example of how even very weak perturbations 
can render the system to become unstable, even if without these perturbations it was perfectly stable. 



VIII. STOCHASTIC SCHRODINGER EQUATION 

Consider the nonstationary Schrodinger equation 

^ = [-iH(r) + af(r,mt)}^, (76) 

in which we set % = 1, = ip(r,a£,t) is a wave function, H(r) is a Hamiltonian, a is real, f(r,t) is a given real 
function, and £(t) is the white noise. With the velocity field defined by the right-hand side of Eq. (76), the Jacobian 
matrix (19) becomes 

J(r, r', at, t) = \-iH(v) + af(r, t)£(t)] 5(r - r') . (77) 
The eigenproblem for the matrix J(a£,t), whose elements are given by Eq. (77), reads 

iK#„ = 4K,t)^. (78) 

Keeping in mind that a is small, the eigenproblem (78) can be solved by means of perturbation theory. In the zero 
approximation, the eigenvector ip n = [ip n (r)] is a column with respect to the spatial variable r, with ip n (r) given by 
the stationary Schrodinger equation 

H{r)i) n {v) =£„Vn(r) , 
so that the zero-order eigenvalue of the Jacobian matrix is 

4°\aZ,t) = -iE n . 
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The first-order approximation for the eigenvalue of the Jacobian matrix is given by 

J n {a^t)=i>+J{a^t)^ n , (79) 

whish yields 

J„«, t) = -iE n + af n (t)£(t) , (80) 

where 

/„(*) = J <(r)/(r,t)^ n (r) dr. 

Note that if /(r, t) = f(t) does not depend on the spatial variable r, then the form (80) with f n {t) — f(t) is an exact 
eigenvalue of the matrix J(a£, t). The multi-index n, labelling the eigenvalues, can be discrete as well as continuous. 

For the stationary eigenvectors ip n of the Jacobian matrix, the multiplier matrix, by theorem 3, possesses the same 
eigenvectors and its eigenvalues are 

Hn t) = exp | -iE n t + a fn(t'W)dt'^ . (81) 
From here, the stability index (47) is 

a(a,t) = a 2 7 /* f n {t') dt' . (82) 
Jo 

The function f(r,t) in Eq. (76) can always be chosen so that to satisfy the inequality 

lim - f f n {t') dt' >0 . (83) 

Switching off stochastic fields results in the neutral motion, for which 

lim lim a(a, t) = . (84) 

But for infinitesimally weak stochastic perturbations, the motion becomes unstable, with 

lim lim a(a,t) = oo , (85) 

where condition (83) is taken into account. In this way, the system described by the Schrodingcr equation is stochas- 
tically unstable, although for some temporal interval, when <j(a, t) <C 1, it can be treated as almost isolated. 

IX. SKETCH OF GENERAL SITUATION 

In the general case, the stochastic field £(t) — [£i(x,t)] is a column composed of the elements &(x, t) depending 
on space as well as on time. This field has to enter the evolution equations as a multiplicative noise. To consider a 
quasi-isolated system, the stochastic term is included with the factor a, which is assumed to be infinitesimally small. 
For a«l, the Jacobian matrix (19) can be calculated by perturbation theory, which yields an expression of the form 

J(a£,t)~J(0,t)+aj'(t,t) . 
In the representation of a basis {(fi n (t)} of vectors ip n (t) — [ip n i{x, t)], this reads 

J mn (a€,t) ~ J mn (0,t) A l mn (x,t)^(x,t) dx . (86) 

i J 

If ip n (t) are the eigenvectors of the multiplier matrix M(a£,i), then its eigenvalues, by theorem 2, are given by Eq. 
(31). Averaging these eigenvalues over stochastic fields, implied to be Gaussian, gives 
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< (j, n (a£, t) > = /i n (0, i) exp 



K„(i) 



(87) 



where the factor 



M„(0,t) = exp|^ J nn {0,t')dt' 



is the multiplier of an isolated system, without any random perturbations, and 

K n(t)=y~] dx 1 dx 2 / A l nn (xi,t 1 )A 3 nn (x 2 ,t2) < £i(xi,ti)£j(X2,t2) > cfei dt 2 

is caused by stochastic perturbations. Then the stability index (47) is 



<r(a, t) = sup Re 



r* a 2 

j Jnn(0,t') dt' + — K n (t) 



(88) 



where a is real. 

When a — > 0, the stability index 



is defined by the properties of the system without perturbations. The limits a 

sup n Re K n (t) 



a(0,t) = sup Re / J nn (0,t') dt' 



lim 



a(0,t) 



oo . 



(89) 

and t — > oo do not commute if 

(90) 



Then the quasi-isolated system is stochastically unstable. This condition is accomplished for the concrete cases 
considered above. It is, of course, impossible to prove that any given quasi-isolated system is, with probability one, 
stochastically unstable. However, the above consideration suggests, with a high level of probability, that there always 
exists such a noise which renders stochastically unstable any particular system. This thesis is certainly correct for 
those systems which, in the absence of noise, display neutral motion. Then Re J nrl (0,t) = 0, hence <r(0,t) — 0, and 
condition (90) is obviously valid. As is shown in Section 7, condition (90) can be hold true even for systems that 
are stable when there is no noise. Let us also emphasize that if, instead of white noise, we would consider infrared 
noise, then condition (90) would necessarily hold. Really, for a deterministic system at large time, one usually has 
cr(0, t) ~ t, while for infrared noise K n (t) ~ t 2 . This makes condition (90) evidently valid. 

Finally, it is important to note that for stochastic dynamical systems the divergence of averaged trajectories is 
not necessarily exponential but may be of algebraic form [4] , which implies that the norm of the averaged stochastic 
multiplier matrix has the power-law behaviour 



< M(a£,t) > || = aA 



(13 > 0) , 



where A is a constant and t c is the chaotization time defining the crossover between stable and chaotic motion. For 
t <C t c , the motion is stable, while for t 3> t Cl it becomes chaotic. The arising instability corresponds to weak chaos 
since the effective trajectory divergence is only algebraic but not exponential. In this case, the local stability index 
(47) is 



<r(a, t) = In 



aA 



From here it follows that the limits 



lim lim a(a,t) = — oo , 



lim lim a(a,t) = +oo 

a— >0 t— >oo 



do not commute with each other. Therefore such a system is also stochastically unstable. 
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X. CONCLUSIONS 



A convenient characteristic for analyzing the stability of dynamical systems is the local stability index (5) . This can 
be expressed through the multiplier matrix M(t) as 

cr(t)=ln||M(i)|| . 

For deterministic (nonstochastic) dynamical systems, there exists another representation of the stability index through 
the Lyapunov or stability matrix ReJ(i), where J(t) is the Jacobian matrix associated with the considered system. 
The name of the Lyapunov matrix comes from the fact that its eigenvalues arc the local Lyapunov exponents. Then 
the stability index, if condition (37) holds, can be written as 

a(t) = f \\Re J(t')\\ alt' . 
Jo 

This presentation, however, is not valid for stochastic dynamical systems. For the latter, the stability index is to be 
calculated by means of Eq. (45). The form (47) of the stability index, 

<r(a,t) = ln|| < M(a£,t) > || , 

is a handy representation for studying the influence of weak stochastic perturbations. The main physical conclusions 
resulting from the general approach and particular examples are as follows. 

(i) Nonexistence of isolated systems. The fact that no real physical system can be completely isolated, but is always 
subject to uncontrollable random perturbations, is more or less generally accepted [5—11]. The point that the concept 
of an isolated system is logically self-contradictory has also been emphasized [9,12,13]. What is principally new in the 
present paper is the demonstration that isolated systems are stochastically unstable with respect to infinitcsimally 
weak random perturbations. A given physical system can be considered as almost isolated, or quasi-isolated, during 
a finite time interval, but it cannot be treated as such for ever. Sooner or later, a quasi-isolated system looses its 
stability. There are no eternally stable systems in nature. 

(ii) Absence of absolute equilibrium. In the theory of dynamical systems, solutions are termed equilibrium if they 
are either constant in time or periodic or quasiperiodic. However, for a quasi-isolated system, no one of these solutions 
can be absolutely stable for infinitely long time. On a finite temporal interval, a solution can correspond to a stable 
equilibrium, but with increasing time, some kind of noncquilibrium behaviour will certainly appear. For instance, 
big fluctuations, driving the system far from equilibrium, may arise [33,34]. Since statistical systems are a particular 
type of real physical systems, they also have to be considered as quasi-isolated. The absence of absolute equilibrium 
for a statistical system implies that large nonequilibrium fluctuations of mesoscopic scale spontaneously appear in the 
system, being randomly distributed in space and in time [35] . If evolution equations do possess an attractor, this has 
to be a chaotic attractor. 

(iii) Irreversibility of time arrow. As far as completely isolated systems do not exist, but there are only quasi-isolated 
systems, the dynamics of such a system, because of the action of random perturbations, can never be reversed so 
that to exactly return to a particular dynamical state. Since quasi-isolated systems are stochastically unstable, any 
trajectory after sufficiently long time will deviate arbitrarily far from the initial point. All that means the irreversibility 
of time. 
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Appendix. 



Here several formulas, related to the Gaussian white noise, are presented, which have been repeatedly used through- 
out the paper. The stochastic variable £(t), corresponding to this noise, centered at zero, has the properties 

« m » = o , « WW) » = 2 7 *(t - f) , 

«C(ti)^ 2 )...C(Wi)» = 0, 

(2n-l)!l 

« £(tl)£(t 2 ) ■ ■ ■ £(*2n) » = (2 7 )" S & - l ^ 5 ^ - *4) • • • S(t 2n -1 ~ *2n) , 

sym 

where J2 sy m implies the symmetrized sum and (2n — 1)!! = (2n)!/2"n! = 1 • 3 • 5 • • • (2n — 1). As an example, a 
symmetrized sum of fij for n = 2 means /12/34 + /13/24 + /i4/23- The integration of over time gives the Wiener 
variable 




For the latter, one has 

« / 2 tu(t) dtu(t) » = 7 (t 2 - ti) , « w 2n+1 (t) » = , « «; 2 "(t) » = ^( 7 t)" . 
In general, any Gaussian variable G{t) satisfies the equality 

< exp G(t) > = cxp I i <G 2 (t)>| . 

For instance, 

<C cxp{aw(t)} 3> = cxp(a 2 7 i) . 

These formulas are sufficient to understand all calculations, related to the averaging over the white noise, which have 
been made in the paper. 
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